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^ ■ Abstract 

> : 

The exact wave functions that describe scattering of a charged particle by a confined 
magnetic field (Aharonov-Bohm effect) and by a Coulomb field are analyzed. It is well 
known that the usual procedure of finding asymptotic forms of these functions which 
admit a separation into a superposition of an incident plane wave and a circular or 
spherical scattered wave is problematic near the forward direction. It thus appears to 
be impossible to express the conservation of probability by means of an optical theorem 
of the usual kind. Both the total cross section and the forward scattering amplitude 
appear to be infinite. To address these difficulties we find a new representation for 
the asymptotic form of the Aharonov-Bohm wave function that is valid for all angles. 
Rather than try to define a cross section at forward angles, however, we work instead 
with the probability current and find that it is quite well behaved. The same is true 
for Coulomb scattering. We trace the usual difficulties to a nonuniformity of limits. 
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1 Introduction 



Not only is 1999 the fiftieth anniversary of Tokyo Metropolitan University, it also marks forty 
years since the publication of the paper of Aharonov and Bohm QXj in which it was shown that 
there are measurable effects that can be attributed directly to the electromagnetic vector 
potential and only non-locally to the magnetic field itself. Of course in non-Abelian gauge 
theories, which are the subject of this conference, the direct importance of the potential does 
not raise any eyebrows, but it did so with respect to electromagnetism at the time. The 
effect was soon confirmed experimentally by Chambers and subsequently by many others 
0. In connection with anniversaries it may be interesting to note that in 1949 (again, fifty 
years ago) Ehrenberg and Siday published a paper @] in which the ray optics of an electron 
moving near a magnetic field were analyzed. It is possible to view this work as a precursor 
to that of Aharonov and Bohm. 

This paper by no means represents a comprehensive review. The literature of the subject 
is enormous. The reader can use the references we have provided to help trace its detailed 
history. We begin in Sec. 2 with a brief review of the Aharonov-Bohm effect. What is 
involved is the effect of the magnetic field on a charged-particle wave function if the particle 
is forbidden from entering the region where the field strength is not zero. We give a heuristic 
derivation of the influence of the field on the interference pattern of a double-slit experiment. 
Aharonov and Bohm showed that if the magnetic field is restricted to an infinite line and if 
the wave function of the particle is chosen to vanish on that line, then the problem is exactly 
soluble. Almost all work on the subject starts from their solution. The wave function is given 
as an infinite sum of partial waves. The bulk of our presentation consists of an analysis of 
this sum. 

In Sec. 3 we apply the usual method of partial wave analysis to calculate the phase shifts. 
This also leads to a formal derivation of the optical theorem. However, since the phase shifts 
are constant for large partial waves the convergence of the procedure is suspect. In Sec. 4 
we review other attempts to find the asymptotic form of the Aharonov-Bohm function and 
to reconcile the results with the optical theorem. Sec. 5 is devoted to a new presentation of 
the Aharonov-Bohm function. It has an asymptotic form which is valid at all angles. In Sec. 
6 we present graphical results of the probability current distribution from which one usually 
calculates cross sections. The behavior near forward angles is quite finite and closely follows 
the double-slit interference pattern mentioned above. 

In Sec. 7 we show that Coulomb scattering shares many of the same features of Aharonov- 
Bohm scattering. Again we provide graphical results of the probability current as derived 
from the exact Coulomb wave functions. 
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2 The Aharonov-Bohm effect and a simple model 



Consider a charged particle in the presence of a magnetic field B(r) which is described by 
a vector potential A(r) such that B = V x A. If in the absence of the magnetic field the 
system is described by a stationary-state wave function ipo{ r ), then in the presence of the 
field it will be described by the path-dependent multiple-valued wave function 

Q 



ifj(r) = ^o(r) exp 



z- / A(r') 

hC Jr Q 



dr' 



V>o(r)e' 



(1) 



Here q is the charge of the particle and P symbolizes the path followed by the integral. 
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Figure 1: The Aharonov-Bohm double slit experiment. 

If one arranges a double slit experiment, as in Fig. 1, so that the part of the beam that 
goes through slit 1 passes entirely above the magnetic field while the part that goes through 
slit 2 passes entirely below the field then using a ray-optics approximation one expects that 
at the particular position that would be the central point of the interference pattern in the 
absence of the field, the wave function will now be 



V(r) 



i^o(r)[e ix(r ' Pl) + e lx(r ' P2) ] 



o(r)e 



*x(r,-Pi) 



1 + exp ( i 



he 



(2) 



where $ is the flux of the magnetic field through the area bounded by the two paths. Thus 
the interference pattern will change with the magnetic field, even if the electron wave function 
is zero in the region to which the field is confined. This contradicts the classical notion that 
the vector potential is basically unphysical and that a charged particle can feel a force only 
at points where the magnetic field itself is nonzero. 

The double slit geometry is not essential for the Aharonov-Bohm effect. Aharonov and 
Bohm treated an idealization of the situation without the slits that could be solved analyti- 
cally. They restricted the magnetic field to be confined to a line of zero width oriented along 
the z direction. That is they took 



B(r) 



z®6(x)d(y). 



(3) 
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In this case, using plane polar coordinates r and 9 in the x-y plane, one may take the vector 
potential to be 



A(r, 



9 



V0. 



2iir 2tt 

The many-valued wave function in the presence of the field would be simply 



ip(r, 9) = ipo(r, 9) exp 



i- 9 

. nc2-K 



M0)e 



iad 



(4) 



(5) 



We have defined a to measure the strength of the flux. It is convenient for later developments 
to take the greatest integer in a to be N and to write 



(6) 



with -~ < 7 < ~. 

The Aharonov-Bohm scattering problem is to find a single- valued solution of Schrodinger's 
equation that corresponds as closely as possible to an incident plane wave moving to the left 
along the x axis from the right of the magnetic field in Fig. 1 (with the source, screen and 
slits eliminated) subject to the boundary condition that the wave function vanishes wherever 
the field is nonzero. For the field of Eq. (|3|) we must take the wave function to vanish at the 
origin. Ignoring the z coordinate, the Schrodinger equation is 



v-if A( r ; 

he 
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(7) 



which simplifies to 
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V>(r,0) = 0. 



The wave function in Eq. (y), which in this case is ip(r,9) = e -' lkrcos0 e ,ia0 ; i s clearly not 
acceptable despite the fact it solves Eq. (R), since it does not vanish at the origin and is 
not single-valued. A single-valued solution can be found by a partial wave expansion and 
becomes a superposition over all / of solutions Ji- a {kr)e lW that are eigenf unctions of the z 
component of angular momentum. The wave function of Aharonov and Bohm is of this form 
and does reduce to the simple plane wave when a = 0. It is 



V>ab(>, 9) 



oo 

E e 

l=—oo 



,iti 



(9) 



Using a concept he calls "whirling waves," Berry || has given a beautiful derivation of this 
result that reconciles the requirement of single-valuedness with the notion that the effect of 
the magnetic field can be given by a simple multiplicative phase factor. 
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3 Two-dimensional partial-wave expansion 



The asymptotic form for large r of a wave function such as that given in Eq. (||) is expected 
to have the form 



-ikr cos 9 



Akr 



+ e i y(8) — 



(10) 



where the scattering amplitude f(9) would have the partial- wave expansion 
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An examination of Eq. (^) shows the phase shifts to be Si = |(i — \l — a\) so that 
2" 



/(*) 




7T 

, sin —a 

irk 2 



N 



l=N+l 



l e iW e -2el _ e -if a 



£ (-i)W 6 ' 

i=— 00 



(12) 



where convergence factors involving the small positive quantity e have been inserted to insure 
convergence. The sums are easily done and the result is 



/(*) 



1 



V2nk 
Letting e — > we find 

/(*) = 



sin 7ra cos | + (1 — cos 7ra)e sin | 
cos 2 § + e 2 



(13) 



sin7ra. <Njv i(N +i\a 1 . 1 — cos7ra r . flN 

-== -1 V (iV+ 5) 9 j. + Z7T (5 COS f , 



(14) 



which is a result obtained by Ruijsenaars 0. Criticism of this formal method of summing the 
partial- wave expansion has been given by Hagen [0. Arai and Minakata || have analyzed 
this expression for f(6) in connection with the optical theorem in great detail without a 
completely satisfactory conclusion. On the other hand Sakoda and Omote have shown that 
the S-matrix is indeed unitary j|. 



4 Integral representations of the Aharonov-Bohm wave 
function 



Another method that makes sense of the sum in Eq. (|9]) is to use a modified version of 
Sommerfeld's representation of the Bessel function 



J v (z) = —e w * I dt e 

Z7T 



-iz cos t ivt 



(15) 
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Figure 2: The contour C of Eq. (|15|). 



where C goes from — tt + rj + ioo to 7r + 77 + ioo with rj a positive infinitesimal as illustrated 
in Fig. 2. 

When this representation is used in Eq. (^) one can perform the sum over / under the integral 
sign to obtain 



4tt 



J_ e {N+\)e 



dt e 



-ikr cos t 



-iyt 



+ 



sin|(t — 6) sin 5^ 



(16) 



One way to proceed is to note that the substitution t — > — t in the second term of the 
integral in Eq. (|l^) leads to an expression identical with the first term but integrated now 
along a contour that is C reflected in the origin. The contours may be distorted, as shown 
by comparison of the left- and right-hand illustrations of Fig. 3, into two straight-line paths 
and a clockwise contour around the pole at t — 9. 



Figure 3: Integration contours for Eq. ( jig 



ikr cos O^iaB 



The contour around the pole gives rise to a modified incident wave of the form e 
which is not single- valued. The other terms may be analyzed by steepest-descent methods 
and give rise to a scattered wave that is valid except for 9 near ±7r. The result is 



i>XB{r,9) 



—ikr cos 9 ^iad — 



1 k r 



-1 



sin an- 



cos \9 



;i7) 



Although the scattered wave is not single-valued, the entire wave function is. It remains 
unclear from this approach, however, how to handle scattering near the forward direction. 
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5 A new integral representation 



We go back to the original contour C and write 



which, by comparison with Eq. fll 



An 



(18) 



defines A(kr,6). By differentiating A(kr,9) with respect 
to kr we can eliminate the singularities. We find, after some manipulation of the resulting 
Bessel functions 

d 



d(kr) 



A(kr,9) = — 2m cos 77? e 



ikr cos 9 



X 



e^V^ftfW (kr) + e-^ e e l ^-^H^_ (kr) 

2+7 2 7 



(19) 



in which Hankel functions appear. Since the wave function must vanish for r = we integrate 
this result from to kr and obtain the following exact expression for the Aharonov-Bohm 
wave function. 
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(20) 



It should be noted that this expression is valid only for — 7r < 9 < 7T. For other values of # 
one should replace 9 by [(6 1 + 7r)mod(27r)] — 71. 

The asymptotic form of ^ab(^, 6 1 ) for large r and any value of 9 may be found by writing 
the integrals from to kr in Eq. (ED) as the difference of integrals from to 00 and from kr 
to 00. For large kr we may use the large-argument asymptotic form of the Hankel functions 
to approximate the second integral (no matter what the value of 9) while the first integral 
can be done exactly |TIJ. The result is 



3 i(jv+i)e -ikr cos e (j~f9 



COS77T 



l-erf(e-^V^cosi0))]}. (21) 



Stelitano [IT] has obtained a similar expression by differentiating the partial-wave sum 
directly with respect to kr. But he treats the / = term separately which means that his 



asymptotic expression is not as compact as Eq. ( 2l[) nor is it strictly correct for backward 



scattering angles. Alvarez |L2[ has also given an asymptotic form for large r as an infinite 
series that is applicable for any 9. 

For y/kr cos ^9 3> 1, (non-forward angles), Eq. fl2~l|) reduces to Eq. (|17|), in agreement 
with the other integration method, while for 9 = n and large kr we get ip(r, n) — > cos cot e 
which is perfectly finite and which exhibits the Aharonov-Bohm interference. 



ikr 
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In view of these results it is to be stressed once again, as has been done by many previous 
workers on this problem, that near the forward direction the total wave function does not 
admit a separation into incident and scattered parts. Hence the usual statement of the 
optical theorem as an expression of unitarity just doesn't make sense. 

6 Properties of the probability current 

From the time-independent point of view the usual optical theorem reflects the conservation 
of the probability current, the integrated form of which is J^ n d9 f • j(r) = 0, for any fixed r. 
The probability current density is given by 

j(r) = -In#*(r)W>(r)]. (22) 
m 

It usually has three terms — coming from the incident wave, the scattered wave and the 
interference between them. In a wave-packet treatment the latter has significance only near 
the forward direction. Because it is proportional to 1/r the scattered term will be quite 
small compared to the incident term. And even at forward angles it will be overwhelmed by 
the incident term. 

The large kr approximation given by Eq. ( plf ) greatly simplifies the numerical calcula- 
tion of r • j(r). The results are plotted in Fig. 4 for 7 = .25 and kr = 100, while in Fig. 
5 we show the naively defined cross section for the same parameters. What we see in the 
Aharonov-Bohm case, where we cannot make the separation into three terms, is that effec- 
tively the forward scattering becomes significant enough to completely modify the behavior 
in the forward direction, while the probability current remains perfectly well-behaved and 
continuous there. 

It is clear from this analysis that the limits of 9 — > 71 and r — > 00 do not commute. It is 
our contention that the former limit is what makes sense physically for any finitely contained 
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Figure 4: A plot of f • j(r) against 9/tt for kr = 100 and 7 = .25, normalized to — 1 at 9 — 0. 
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Figure 5: A plot of the usual differential cross section in units of 1 /k against 9/rr for kr = 100 
and 7 = .25. 



experiment. 



7 Coulomb scattering 

A similar nonuniformity of limits appears in the consideration of Coulomb scattering. The 
exact wave function, written in terms of a confluent hypergeometric function is 

^coui(r) = e lkrcose ^(-irn l;iAx[l - cosfl]), (23) 

with Coulomb parameter 77 = mZe 2 / (hk). For large kr(l — cos#) this has the asymptotic 
form 

Jkrcose e ikr e~ 2 W r(l + «7). 



^coui(r) ~ e i,?logfer(1 " cose) 



" kr (I- cos 9y+ 2i nr (I -iri) V 



(24) 



which looks like a plane wave incident along the z axis plus a spherical wave, but with 
Coulomb distortion. The "r-dependent" scattering amplitude extracted from this is 

V T(l + iri) e-^\ozkr 

JK ' kvii-irj) (i-cos^) 1+2 ^ 1 ' 

and the differential cross section is 

da rf 



dVl k 2 {l -cos9) 2 ' 
When extrapolated to forward angles these are both infinite. 



(26) 



On the other hand, the exact wave function and the probability current are quite finite 
at — 0. In fact the probability current becomes very small compared to its maximum value 
at near-forward angles. This is illustrated in Fig. 6 with r\ = 1.5. 
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Figure 6: Radial component of the probability current for Coulomb scattering, plotted 
against 8/ir and normalized to —1 at backward angles. (The forward angle is now at 9 = 0). 
The parameters are kr = 100 and rj = 1.5. 
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